THE DETECTION OF COLD DUST IN CAS A: EVIDENCE FOR THE 
FORMATION OF METALLIC NEEDLES IN THE EJECTA 
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ABSTRACT 



■ Ejecta from core collapse supernovae contain a few solar masses of refractory el- 

ements and can therefore be the most important source of interstellar dust if these 
elements condense efficiently into solids. However, infrared observations of young su- 
, pernova remnants, such as Cas A or Kepler, and observations of SN 1987A have detected 

'sj" I only ~ 10~^ M0 of hot dust in these objects. Recently, Dunne et al. (2003) obtained 

O 



450 and 850 /um SCUBA images of Cas A, and reported the detection of 2-4 Mq of cold, 
18 K, dust in the remnant. Here we show that their interpretation of the observations 
faces serious difficulties. Their inferred dust mass ignores the effect of grain destruction 



o 

C3 ■ by sputtering, and is larger than the mass of refractory material in the ejecta of a 10 to 

30 M.Q star. The cold dust model faces even more difficulties if the 170 //m observations 
of the remnant are included in the analysis, decreasing the cold dust temperature to 



O I ~ 8 K, and increasing its mass to > 20 M©. We offer here a more plausible interpretation 

I of their observation, in which the cold dust emission is generated by conducting needles 

I in the ejecta. The needle properties are completely determined by the combined sub- 

millimeter and X-ray observations of the remnant. The needles are collisionally heated 
^ ' by the shocked gas. They are very efficient emitters at submillimeter wavelengths, and 

^ . with a resistivity of a few cm can readily attain a temperature of 8 K. Taking the 

destruction of needles into account, a dust mass of only 10~^ to 10^'^ Mq is needed 
to account for the observed SCUBA emission. The needles consist of metallic whiskers 
with < 1% of embedded impurities, that may have condensed out of blobs of material 
that were expelled at high velocities from the inner metal-rich layers of the star in 
an asymmetric explosion. Conductive needles may also be the source of the cold dust 
emission detected by Morgan et al. (2003) in Kepler. Aligned in the magnetic field, 
needles may give rise to observable polarized emission. The detection of submillimeter 
polarization will therefore offer definitive proof for a needle origin for the cold dust 
emission. Supernovae may still be proven to be important sources of interstellar dust, 
but the evidence is still inconclusive. 



Subject headings: ISM: supernova remnants - ISM: individual (Cassiopeia A) - 
ISM: interstellar dust - Infrared: general 



- 2 - 



1. INTRODUCTION 

Cas A, the young remnant of a Type II supernova (SN), has been suggested as an ideal 
observing target for dust that may have formed in the explosive ejecta (Dwek & Werner 1981). 
Initial ground-based searches for dust (Wright et al. 1980, Dinerstein et al. 1982) yielded negative 
results. Dust was first detected in Cas A in the Infrared AstronomAcal Satellite (IRAS) all sky 
survey (Mezger et al. 1986, Dwek et al. 1987, Braun 1987, Arendt 1989, Sakcn, Fesen, & Shull 
1992). However, due to the low spatial resolution of the ZR^S" images, its origin whether interstellar, 
circumstellar or supernova, could not be ascertained. Higher resolution infrared images of Cas A, 
obtained with the Infrared Space Observatory (ISO) satellite, revealed the presence of hot dust in 
the fast-moving knots of the remnant, thereby unambiguously establishing its SN origin (Lagage 
et al. 1996, Arendt, Dwek, & Moseley 1999, Tuffs et al. 1999). Dust temperatures were typically 
~ 100-200 K, consistent with those expected for shock-heated dust. Dust masses were between 
10^^ — 10~^ Mq, depending on the size of the observed region, significantly below the mass of 
condensable elements expected to be present in the ejecta of a massive progenitor star (Woosley &: 
Weaver 1995). 

Supernovae are potentially the most important source of interstellar dust and they play an 
important role in establishing the elemental depletion pattern in the interstellar medium (Dwek 
& Scalo 1980, Dwek 1998, Jones 2000). Furthermore, the early rise in the UV-optical opacity of 
high rcdshift galaxies suggests copious dust production in these objects (Todini & Ferrara 2001). 
The lack of any observational evidence for the formation of massive amounts of dust in supernovae 
is therefore quite frustrating. In light of this, the recent detections of emission from cold dust in 
Cas A (Dunne et al. 2003) and Kepler (Morgan et al. 2003) are very encouraging, since they may 
provide the first observational confirmation of the role of supernovae as prominent interstellar dust 
sources. Their findings require therefore careful scrutiny. In our analysis we will concentrate on 
the cold dust detection from Cas A. 

The 450 and 850 fim SCUBA images of Cas A (Dunne et al. 2003) revealed a new, cold, dust 
component, primarily concentrated in the southern half of the remnant, and bounded between the 
forward and reverse shocks, located at angular radii of, respectively, 153 "and 95 "(Gotthelf et al. 
2001). The total IR spectrum from Cas A, comprised of the IRAS (12 - 100 ^um) and SCUBA 
observations, from which the synchrotron component has been subtracted, was fitted by Dunne et 
al. (2003) with a hot (112]']2i ^) ^ cold (18l4g dust component. The hot dust component 
detected by the IRAS has a mass of ~ 10~^ Mq, and consists of supernova dust heated by the 
interaction of the ejecta with the ambient medium (Arendt, Dwek, Sz Moseley 1999). The cold dust 
has a mass between 2 and 12 Mq, depending on its radiative properties (Dunne et al. 2003 and 
§2 below). 

The properties and origin of the cold dust are difficult to explain even if one adopts the lower 
mass estimate. The inferred dust mass is too large to be of either circumstellar or interstellar origin. 
Such origin will suggest the presence of an associated ~ 200 Mq mass of hydrogen, clearly ruling 
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out the former. An interstellar cloud of that mass interacting with the remnant will give rise to 
optical and IR emission lines and to dynamical effects that are unobserved. Below we argue that 
the inferred dust mass is also too large to be of supernova origin. Furthermore, the gas associated 
with the dust, will collisionaly heat the dust to a temperature that is significantly higher than that 
inferred from the SCUBA observations. Nevertheless, ruling out an interstellar or circumstellar 
origin, Dunne et al. (2003) adopted the least implausible explanation, and concluded that the dust 
must be of SN origin. 

This interpretation of the SCUBA observations faces several difficulties: 

1. The required dust condensation efficiency in the SN ejccta is too high even for the amorphous 
or clumpy dust aggregates. Cas A is believed to be the remnant of the explosion of a 10 to 
30 M0 star (Fesen et al. 2001). The total mass of condensible elements (C, Mg, Si, S, Ca, Fe, 
Ni, and associated O) ejected in the explosion of 10 and 30 M© stars is about 0.3 and 2 Mq, 
respectively (Woosley & Weaver 1995). So even the lowest estimate of ~ 2 Mq of cold dust 
requires a very unlikely condensation efficiency of 100% in the supernova ejecta of a 30 M0 
star. The large dust mass becomes even more embarrassing if, as suggested by the distribution 
of elemental abundance ratios in the ejecta, the progenitor of Cas A was only a 12 Mq star 
(Willingale et al. 2002). Even allowing for a massive progenitor, the required condensation 
efficiency is too high considering the fact that important refractory elements (Si, Ca, and Fe) 
are observed to be in the gas phase of the remnant, and to be more uniformly distributed 
than required by the Dunne et al. (2003) interpretation of the SCUBA data (Hwang, Holt, 
& Pctrc 2000; Hughes et al. 2000, Willingale et al. 2002). Finally, the total mass of dust 
that must have formed in the ejecta may actually be greater than implied from the SCUBA 
observations, since a fraction of it may have been destroyed by thermal sputtering. 

2. The density of the gas associated with the cold dust is too high to account for the low dust 
temperature. Associated with the refractory elements produced in a 30 Mq star are about 
4 Mq of other abundant metals, primarily oxygen, that are not locked up in dust. We 
can derive a lower limit to the electron density of the emitting region by assuming that it 
is smoothly distributed between the forward and reverse shocks, and that it covers about 
15% of the sky as seen from the center of the explosion. The volume of the region, for an 
adopted distance of 3.4 kpc, is then ~ 2 x 10^^ cm~^, giving a minimum electron density of 
~ 3 cm~^. Actual electron densities are significantly higher, since oxygen is expected to be 
multiply ionized in the emitting region. However, even at the low density of 3 cm^^, typical 
0.01 /xm-sized supernova condensates will be coUisionally heated to temperatures of about 
50 K (Dwek 1987). Put differently, a density of 3 cm~^ is about 100 times higher than that 
required to heat the dust particles to the observed temperature of 18 K. To achieve this cold 
dust temperatures at these high electron densities, the radii of the dust particles in the ejecta 
must be > 1 ^m. However, such radii are much larger than those predicted from nucleation 
calculations (Kozasa, Hasegawa, & Nomoto 1991). 
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3. Finally, the 170 /xm ISO observations of Cas A are actually a detection, implying an even 
larger amount of colder dust. In their analysis, Dunne et al. (2003) regarded the ISO 170 ^m. 
observations of Cas A (Tuffs ct al. 1999) as a lower limit, predicting a significantly higher 
flux from the remnant than observed at that wavelength. However, it was pointed out by the 
referee that it is difficult to explain how the ISOPHOT instrument on board the ISO could 
have missed most of the emission from that region of the remnant. In fact, the observed 
170 ;um emission is an actual detection (private communications between the referee and R. 
J. Tuffs made available to the author), and the 170 ^m. flux can be entirely attributed to 
the hot dust emission. If the 170 /ini flux produced by the cold dust is constrained to be 
less than 10 Jy, then the SCUBA data can only be fitted with ~ 7 to 9 K dust assuming a 
emissivity law. The mass of cold dust is then increased to ~ 24 to 16 Mq, respectively, 
even if it has the high submillimeter radiative efficiency of amorphous or clumpy aggregate 
material. So regarding the 170 /im flux as an upper limit on the emission from the cold dust 
component simply exacerbates the difficulties listed above. 

The problems listed above can be alleviated if the dust responsible for the observed emission 
consists of conducting needles or whiskers. The high conductivity (low resistivity) of these needles 
gives rise to a mass absorption efficiency that is several orders of magnitudes greater than that of 

interstellar graphite or silicates, and even the amorphous clumpy aggregates considered by Dunne 
et al. (2003). Sufficiently large needles are efficient antenna-like radiators out to millimeter wave- 
length, depending on their length-to-diameter ratio. Their large radiative efficiency significantly 
lowers the mass of cold dust required to explain the SCUBA observations. It also increases the 
density of electrons needed to collisionally heat the dust to the observed temperature. The intro- 
duction of such efficient radiators as the source of the observed submillimeter dust emission requires 
therefore a balance between the need to reduce the mass of the emitting dust without increasing 
the plasma density beyond the limits that may be imposed from X-ray observations. 

Conducting needles have been considered before, in a different context, by Hoyle &: Wick- 
ramasinghe (1999) and Wickramasinghc (1992) as a source of starlight opacity that could have 
produced a non-cosmological cosmic microwave background (CMB) in a steady state cosmology. 
They were also considered as a way to thermalize starlight and generate the CMB after the big 
bang in cold big bang models (Wright 1982, Li 2003). It should not be too surprising to find needles 
in supernova ejecta. Supernovae condensates may originally be different in morphology or compo- 
sition from the dust particles commonly found in the diffuse interstellar medium (ISM), and arc 
likely to undergo additional chemical and physical processing in the ISM. There is prior evidence 
that SN condensates may be different from their ISM counterparts. The hot dust emission from 
Cas A detected by the ISO was interpreted as emission from protosilicate material (Arendt, Dwek, 
k. Mosleley 1999), and conducting needles may have been detected in SN 1987a (Wickramasinghe 
&; Wickramasinghe 1993). 

In the model proposed here, the metallic needles were formed in metal-rich ejecta that were 
expelled at high velocities into the ambient medium. The needles are collisionaly heated by the 
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reverse shock resulting from the interaction of the ejecta with the ambient medium. In the fohowing 
section (§2) we examine the different dust component contributing to the observed IR to submil- 
limetcr spectrum of Cas A, showing that conducting needles arc a necessary emission component if 
one wants to fit the 170 //m ISO observations with a reasonable dust mass. In §3 we calculate the 
properties of the needles and the electron density required to explain the SCUBA observations for 
various type of conducting material. The results of the paper are summarized and discussed in §4. 



2. DUST EMISSION COMPONENTS IN THE SPECTRUM OF CAS A 

Figure 1 depicts the observed infrared (IR) to submillimeter spectrum of Cas A comprising of 

the 12, 25, 60, and 100 /um data obtained by the Infrared Astronomical Satellite (IRAS) (Dwck ct 
al. 1987), the 170 fim datum point obtained by the Infrared Satellite Observatory (ISO) satellite 
(Tuffs et al. 1999), and the 450 and 850 fim SCUBA data (Dunne et al. 2003). An almost identical 
850 iim flux was reported by Loinard et al. (2003). Also shown in the figure is the synchrotron 
emission spectrum, extrapolated with a power law (Mezger et al. 1986). Table 1 summarizes 

the data used in this paper. 

The specific flux at wavelength A received from an ensemble of dust particles with a total mass 
Md, located at a distance D, and radiating at an equilibrium dust temperature T^, is given by: 

FuW = (^) '^(A)i?.(A, T,) (1) 

where k{X) is the mass absorption coefficient of the dust, and Bi, is the Planck function. For a 
k{\) oc A"" emissivity law, eq. (1) can be written as: 



F^{X) = 8.30 X 10 
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^kpc 



cw? g ^ 



A-3 

y 1 



^ / 14387.7 _ ^ 

\ XnroTd 



-1 



Jy (2) 



where k{\q) is the mass absorption coefficient at some reference wavelength Aq 



Regarding the ISO observations as a lower limit, Dunne et al. (2003) fitted the Cas A spectrum 
with a hot and cold dust component radiating at an equilibrium temperature of 112 K, and 18 K, 
respectively. Figure 1 depicts their 2-component fit to the IRAS and SCUBA data of Cas A. 
Interstellar dust particles have a value of «;(450 /um) ~ 2.5 cm^ g~^ (Draine & Lee 1984). Assuming 
that the Cas A dust has similar radiative properties gives a total mass of ~ 2 X 10~3 Mq of hot 
dust and about 10 of cold dust. The mass of cold dust may be lower if it has an amorphous 
or clumpy aggregate structure, for which laboratory measurements and theoretical models give a 
mass absorption coefficient of k(450 jim) = 15 cm^ g~^ (Ossenkopf & Henning 1994). The dust 
mass is then only about 2 M©, still incompatible with a circumstellar or interstellar origin, but 
more likely (according to Dunne et al. 2003) to be of supernova origin. 

The model of Dunne et al. faces more difficulties if the ISO 170 /xm observations are an actual 
measurement (Tuffs et al. 1999). The observed flux can then be entirely attributed to the emission 
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Fig. 1. — The infrared-submillimeter spectrum of Cas A consisting of the IRAS 12, 25, 60, and 
100 /im (filled triangles, Dwek et al. 1987), the ISO the 170 fxm (filled diamond, Tuffs et al. 1999), 
and the SCUBA 450 and 850 /im observations (open squares, Dunne et al. 2003). Regarding the 
170 /im fiux as a lower limit, Dunne et al. fitted the data with a hot+cold dust components. 
Taking the 170 jim flux as a detection, the data can be fitted with a hot dust+conductive needles 
components. Needles are required to alleviate the excess mass problem. 
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from the hot dust component (see Figure 1). A fit to the SCUBA observations would then lead to 
a colder dust temperature and therefore to a significantly higher dust mass, even if it consists of 
amorphous or clumpy aggregates. 

The mass of cold dust depends on the mass absorption coefficient, and can be further reduced 
if the dust is a more efficient emitter at far-IR and submillimeter wavelengths than the amorphous 
clumpy aggregates considered by Dunne et al. (2003). Conducting needles are just such kind of 
dust. 



3. CONDUCTING NEEDLES 

3.1. Mass Absorption Coefficient 

Consider a conducting needle-like dust particle represented by a circular cylindrical of radius 
a, length i (a <^ i), mass density pm, and mass = TTO^ipm- Most of the incident radiant energy 
is absorbed by ohmic dissipation. The absorption cross section is given by aabs = P/S, where 
S = (c/87r)£'^ the time-averaged Poynting vector of the incident radiation, E the amplitude of the 
incident electric field, and P the power absorbed in the grain, are given by (Wright 1982): 

^ ~ 2R~3 2R 

l7ra%E£)2 

3 2pr£ 

In the equation above, V is the voltage on the needle, R = pui/TTo^ is its resistance, p^ its resistivity, 
and the factor of ^ arises from averaging over the angles of incident between the electric field and 
the needle (Wright 1982). The absorption cross section is aabs = (47r/3c)7ra^£/pi^, and the mass 
absorption coefficient is: 

1^0 = = -IT ■ (4) 

ma 3 \ PmPRC J 

The mass absorption coefficient of the needles is nearly constant for wavelengths > 100 /xm 
(Hoyle &; Wickramasinghe 1999, Wright 1987), with a long-wavelength cutoff at Aq given by (Wright 
1982): 

after which k(A) oc at wavelengths > Aq. 

Table 2 lists several physical properties of different conducting and semi-conducting material 
consisting of elements that are likely to condense out of metal-rich SN ejecta. The electrical 
resistivity of metals is generally a sum of two independent contributions: (1) the resistivity caused by 
the collisions of the conduction electrons with lattice phonons, and (2) the resistivity caused by the 
scattering of the electrons by impurity atoms and lattice defects. The first contribution is strongly 
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temperature dependent, generally dropping by 2 to 3 orders of magnitude from room temperature 

to 20 K. For example, the resistivity of pure iron drops from ~ 10~^ J7 cm at 293 K to a value of 
~ 3 X 10~® r2 cm at 20 K (1 $7 cm = 1.139x10^^^ sec). The contribution of impurities and lattice 
defects to the resistivity, also referred to as the residual resistivity of the material, is independent 
of temperature, and can dominate the resistivity at sufficiently low temperatures. For alloys, such 
as nichrome (80% Ni-20%Cr), the resistivity is virtually constant between temperatures of 300 and 
20 K. At very low concentrations, impurities can contribute about 4 /iQ cm per percent of impurity 
to the resistivity (Kittcl 1963). SN condensates arc likely to have embedded impurities, however, 
their concentration is unknown. Given the uncertainties in the resistivities of some elements at 
very low temperatures, we will take the value of pn as a variable in our calculation. 



3.2. Temperature and Long Wavelength Cutoff 

For given material properties, the needle emission spectrum is completely determined by its 
temperature and long- wavelength cutoff. The observational constraints are the Fi,(450 pm) to 
Fiy{850 nm) flux ratio, which is equal to =3.5ib2, and the rapid drop in the spectrum at 170 /xm. 
Figure 2 depicts the 450-to-850 fim (dashed lines) and the 450-to-170 ^m flux ratios (solid lines) 
as a function of the cutoff wavelength and needle temperature. The figure shows that the SCUBA 
observations arc well matched by 8.2 K needles with a long wavelength cutoff of Ao = 400 fim. 
Figure 3 depicts the £/a ratio for Aq = 400 /Ltm as a function of needle resistivity. The £/a ratio 
varies from ~ 4000 to about 100 over the range of resistivities depicted in the figure. For nichrome, 
which has a constant resistivity as a function of temperature, the £/a ratio is about 260 giving 
a needle length of ~ 2.6 /um for a radius of a = 0.01 /um. Figure 1 depicts the fit of the needle 
spectrum to the SCUBA data. The needle model predicts a significantly lower flux at wavelengths 
below ~ 300 nm compared to the cold dust model of Dunne et al. (2003). The 100 to 400 /xm flux 
from Gas A is therefore an important discriminator between the cold dust and the needle model. 



3.3. Total Mass 

The total mass Ma, of needles radiating at temperature is given by: 

where S^{X) is the observed flux density, and 15=3.4 kpc is the distance to the remnant. Figure 3 
shows the mass of dust required to produce a 450 /um flux of 45.4 Jy for a needle temperature of 
8.2 K. An average material density of 8.15 g cm~^ was adopted in these calculations. Unlike the 
cold dust model, the mass of elements locked up in needles is only a small fraction of the total mass 
of metals produced in the explosion. However, the masses shown in the figure are only lower limits 
on the mass of dust that condensed out of the ejecta, since the needles are destroyed by thermal 
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Fig. 2. — Contour levels of the Ri = Fjy(450 fim)/Fi,{850 /iin) (dashed lines) and the R2 = 
Fi/(450 /nm)/F;y(170 /xm) (solid lines) flux ratios are plotted versus the cutoff wavelength, Aq, and 
needle temperature. A good fit to the observational constraints Ri = 3.5 it 2 and R2 » I require 
Ao = 400 fim, and a needle temperature of 8.2 K. The needle spectrum is depicted in Figure 1. 
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sputtering in the shock-heated gas. In the following we calculate the density of the gas required to 
heat the needles to a temperature of ~ 8 K, and the corresponding rate of grain destruction. 



3.4. The Heating and Cooling of Needles 

The heating rate of a needle due to electronic collisions is given by: 

HdiT) = 27ra£ne{l-fr) J 9{E)v{E)Edep{E)dE 



(7) 



= 2mdne(l - fr) J g{E)v{E) 



EdepjE) 

. Pma 



dE 



where g{E) = 2'K~^/'^{kT)~^/'^E^/'^ eyq){—E /kT) is the Maxwell distribution function of electron 
energies E at temperature T, v{E) and rig are, respectively, the electron velocity and number 
density, is the fraction of incident electrons that are reflected from the grain surface, and E^epiE) 
is the energy deposited by incident electrons with energy E in the needle. If the electrons are stopped 
in the grain, E^ep = E, otherwise E^ep = x{dE/dx), where x pa 2a is an average distance traveled 
by an electron through the needle, and (d£^/da:;) is the electronic stopping power (the energy lost 
in the solid per unit length traversed by the electrons). 

The range of electrons with energies between 370 and 10^ eV is given by (Iskef, Cunningham, 
& Watt 1983): 



i?e(A*m) = 2.59 X 10" 



g cm 



E{eV) 



1.492 



(8) 



Taking R^, = 2a = 0.02 //m as the effective pathlength through the cylinder, we get that the 
maximum energy at which the electrons are stopped in the needle is given by: 



E^{eV) = 403 



-[ 0.67 



1.64 X 10' 



cm 

,3 



(9) 



for pr, 



.15 gr cm 



Differentiating eq. (8) w.r.t. E gives the electronic stopping power in the 370 to 10^ eV energy 
range: 



ax \g cm \e\ J 

The dust heating rate can be now written as the sum: 



-0.492 



eV cm ^ 



(10) 



Hd{T) = 2mdne{l- fr){ r g{E)v{E)(^\dE (11) 

Jo \PmO-J 
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Radiating at temperature Tfi, the needle cools at a rate given by: 

Ld = 4:md J n{vQ)TTB^{u, Ta) diy (12) 

where 

k{X) = kq for A < Ao (13) 

-2 

for A > Ao (14) 

where kq and Aq are given by eqs. (4) and (5), respectively. 

In equilibrium, the dust cooling rate equals its heating rate by electronic collisions. Figure 4 
depicts the electron density ng as a function of needle resistivity at various fixed plasma tempera- 
tures, represented by the diagonal lines in the figure. The figure depicts two plasma components 
identified in the X-ray analysis of Cas A (Willingale et al. 2003) : a hot component characterized by 
{ue, kTe} = {16±3 cm-^, 3.27±0.86 keV}, and a cool component with {ue, kTe} = {61±15 cm"^, 
0.45ib0.10 keV}. The shaded areas in the figure represent the allowable density and temperature 
combinations for the two plasma components, defined by the la uncertainties around their central 
values. In the hot plasma, the needles are mostly transparent to the incident electrons. The dust 
heating rate is then independent of electron temperature (Dwek 1987), accounting for the narrow- 
ness of the region between the two temperatures that bound the hot plasma. The figure shows that 
the needle resistivity needs to be between ~ (2 — 5) x 10~^ Q cm, and could reside in either the 
hot or cool (or both) plasma components of the remnant. 




3.5. Needle Destruction by Sputtering 

Dust particles in a shocked plasma are eroded by thermal sputtering. For plasma temperatures 
of ~ 10^ K (~ 0.86 keV), the lifetime of silicate dust against sputtering by heavy nuclei, primarily 
oxygen, in the SN ejecta is given by (Dwek, Foster, &: Vancura 1996): 

r(yr)»5x 10=^1=1 (15) 

where n{0) is the number density (in cm~^) of the oxygen nuclei in the hot gas. Iron dust has a lower 

surface binding energy of about 4.3 cV, compared to the surface binding energy of silicates, which 
is about 6 eV (Dwek & Scalo 1980). Iron dust will therefore have a shorter lifetime than silicate 
particles of identical size. To be definite, we adopt a value of 10"'' for the numerical coefficient in 
the expression for the dust lifetime. Assuming that the density of heavy elements is approximately 
given by the electron density in the two different plasma components, and taking a = 0.01 fxm, 
we get that the lifetime of iron needles is about 60 yr in the cool, and about 16 yr in the hot 
plasma component. Both times are shorter than the range of ages of these two plasma components 
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Fig. 3. — Plot of the mass of 8.2 K needles required to fit the SCUBA data (see Figure 1), and the 
needle length-to-radius, £/a, ratio as a function of needle resistivity. Also shown in the figure are 
the room temperature resistivities of iron (Fe), and nichrome (Ni-Cr). 
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Fig. 4. — The electron density required to heat a needle to 8.2 K is plotted versus needle resistivity 
for different plasma temperatures. The shaded areas define the parameters of the cool and hot 
plasma components in Cas A, defined by {ng, Tg} = {61±15 cm~^, 0.45ib0.10 keV}, and {16ib3 
cm~^, 3.27±0.86 keV}, respectively. 
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which are, respectively, about 80 (20-273) and 131 (101-182) yr (Willingale et al. 2003). The 
radiating dust mass represents therefore the fraction of the total mass of shocked dust which dust 
has not yet been destroyed. A rough estimate for the actual mass of needles that must have been 
condensed in the ejecta in order to produce the observed emission can be obtained by multiplying 
the masses given in Figure 3 by the ratio of the age of the plasma to the sputtering lifetime. Typical 
needle masses need therefore to be 2 to 10 times larger than indicated in the figure if they reside, 
respectively, in the cool or hot plasma component of the remnant. Adopting average resistivities 
for the needles in the two plasma components (Figure 4), we get that the total mass of needles 
required to produce the SCUBA observations is ^ xlO~^ Mq if they reside in the cool plasma, and 
~ X 10~^ M0 if they reside in the hot plasma. There may yet be an additional reservoir of colder 
dust. So even after the effects of grain destruction have been taken into account, the total mass of 
needles required to have formed in the ejecta still represents a very small fraction of the mass of 
available metals in the ejecta. 



4. SUMMARY AND DISCUSSION 

Recent SCUBA observations of Cas A (Dunne et al. 2003, Loinard et al. 2003) have revealed 
an excess of submillimeter emission above the extrapolated synchrotron emission from the remnant. 
Dunne ct al. (2003) interpreted the excess as emission from about 2 Mq of cold, 18 K, dust that 
has formed in the SN ejecta. This dust mass requires the implausible scenario that essentially all 
the condensible elements that formed in ejecta be locked up in the dust. X-ray observations rule 
out such scenario. Their interpretation faces more difficulties if grain destruction by sputtering is 
taken into account, and if the 170 //m observations of the remnant are also included as a detection 
in the analysis. The cold dust temperature then becomes 8 K, increasing the inferred cold dust 
mass to about 20 Mq. 

A conductive needle origin for the cold dust emission from Cas A provides a simpler alternative 
explanation than the one put forward by Dunne et al. (2003), even when the 170 /im observations 
are included in the analysis. The results of the paper can be briefly summarized as follows: 

1. the submillimeter spectrum can be fit with 8.2 K needles having a long wavelength cutoff of 
400 /xm (Figures 1 and 2), 

2. adopting a range of 10~^ to 10~^ Q, cm of material resistivities, the mass of needles required 
to produce the SCUBA observations falls between 2xl0~^ and 0.02 Mq, respectively, and 
that the needle length-to-radius {£/a) ratio is between 4000 and 100, respectively (Figure 3), 

3. balancing the radiative cooling rate of the needles at 8.2 K to the coUisional heating rate 
for a range of plasma densities and temperatures, the range of needle resistivities that are 
compatible with observational X-ray constraints can be narrowed to (2-5) x 10~^ Q, cm (Figure 
4). These values suggest that the needles are made of pure conducting material (iron?) that 
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includes about 0.5 to 1% of embedded impurities 

4. the needles are destroyed by sputtering in the hot gas. The needle lifetime is shorter than 
the ionization age of the various X-ray emitting components, suggesting a continuous supply 
of newly swept-up dust. The mass of needles required to produce the SCUBA observations 
is therefore increased to ~ 10^'^ Mq if they reside in the cool plasma, and ~ 10~^ Mq if they 
reside in the hot plasma. There may be an additional reservoir of cold dust that has not been 
shocked yet. 

The needle model proposed here is consistent with a scenario in which the metallic needles 
are formed in bullets of metal-rich ejecta that were expelled at high velocities into the ambient 
medium. Such scenario is similar to the one observed in SN 1987A, in which the increased ionization 
state of its envelope, the early appearance of X- and 7-rays from the radioactive decay of ^^Co, 
and the width of the spectral lines of heavy elements, suggested that metal-rich material from the 
inner region of the ejecta was macroscopically mixed and expelled at velocities of ~ 3000 km s~^ 
(Graham 1988, Arnett et al. 1989). The same picture emerges from optical and X-ray studies of 
the ejecta of Cas A. The spatial distribution of the nucleosynthesis products in the remnant shows 
that the most Fe-rich material is very clumpy and located at the outer edge of the ejecta (Hughes 
et al. 2000, Willingale et al. 2002), suggesting that the stellar layers were macroscopically mixed, 
spatially inverted, and ejected as bullets into the ambient medium in an asymmetric explosion. 
The interaction of the ejecta with the ambient medium will drive a reverse shock wave through 
the ejecta. The shocked plasma will then collisionally heat the dust to the observed temperature. 
The submillimeter emission detected by the SCUBA from Kepler (Morgan et al. 2003) may have 
a similar origin. 

Conductive needles may be efficiently aligned in a magnetic field, so that the cold dust emission 
may be partially polarized, depending on the geometry of the magnetic field. The submillimeter 
polarization from the dust can dominate that from the synchrotron emission. The detection of 
polarization at wavelengths < 450//m, and a drop in the cold dust spectrum from 400 to 170 /xm 
will therefore be definitive confirmations of the proposed model. 

The results of our analysis show that the SCUBA observations need not imply the detection of 
massive amounts of dust in the remnant. Conductive needles offer a better explanation with only 
a fraction of the mass inferred by Dunne et al. (2003) and Morgan et al. (2003). Nevertheless, 
supernovae may yet be important sources of interstellar dust. Here we argue that the evidence put 
forward sofar is still inconclusive. 
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Table 1. Observed Infrared and Submillimeter Fluxes Prom Cas A 



Emission component 


12 fim 


25 /xm 


60 fim 


100 /Ltm 


170 nm 


450 /xm 


850 urn 


Total^ 


16.9 ±2.5 


152±23 


123±20 


71±20 


35±10 


69.8±16.1 


50.8±5.6 


Synchrotron^ 


2.0 


3.3 


6.1 


8.7 


12.5 


24.4 


37.8 


Ionic lines 


4.3 


< 2 


< 15 


< 1 


11.5 






Thermal dust 


10.6±2.5 


149±23 


117±20 


62±20 


lliblO 


45.4±16.1 


13.0±5.6 



^Fluxes in units of Jy. IRAS 12, 25, 60, and 100 fim data arc from Dwek ct al. (1987), 170 fiui 
flux from Tuffs ct al. (1999), and SCUBA 450 and 850 ^m fluxes arc from Dunne ct al. (2003). 

^Synchrotron fluxes calculated for a power law (Mezger et al. 1986), normalized to 2.0 Jy 

at 12 jj,m. 



Table 2. Summary of Physical Properties of Different Materials 





Pure iron 
(Fe) 


Nichrome 
(Ni-Cr) 


Silicon carbide 
(SiC) 


Graphite 
(C) 


Mass density, pm(g cm~^) 


7.8 


8.5 


3.2 


2.2 


Resistivity, pr{Q cm) 


9.7x10"^ 


1.1x10"* 


2.0x10"'' 


9.1x10"" 


Mass absorption coefEcient, Ko{cm^ g~^) 


1.3x10'' 


l.OxlO'^ 


1.5x10^ 


4.8x10* 



^Resistivities measured at 293 K. Data taken from the ORG handbook of Chemistry and Physics, 84th Edition, 
2003-2004. For unit conversion: 1 Q cm = 1.139x10"^^ s. The value of ko was calculated using eq. (4). 



